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Abstract 

A new type of cosmic-ray experiment is under construction in the Pyhasalmi mine 
in the underground laboratory of the University of Oulu, Finland. It aims to study 
the composition of cosmic rays at and above the knee region. The experiment, called 
EMMA, will cover approximately 150 m 2 of detector area. The array is capable of 
measuring the multiplicity and the lateral distribution of underground muons, and 
the arrival direction of the air shower. The full-size detector is expected to run by 
the end of 2007. 

Key words: cosmic rays, composition, muon lateral distribution, muon 
multiplicity, air shower, underground experiment 


* Corresponding author. 

Email address: timo.enqvist@oulu.fi (T. Enqvist). 


Preprint submitted to Elsevier Preprint 


5 February 2008 





1 Introduction 


Due to a slight change observed in the cosmic-ray energy spectrum in the 
energy interval of 10 15 - 10 16 eV, it is believed that the origin, modification in 
the chemical composition, acceleration mechanism or propagation of cosmic 
rays (or a combination of these) changes. Up to this energy, so called knee 
region, most cosmic rays are supposed to be produced inside the galaxy, and 
are also confined by the galactic magnetic held. 

At these high energies the source cannot be observed directly and the cosmic- 
ray composition is used as a tool to investigate the origin of the cosmic radia¬ 
tion. The direct composition measurements are no longer practical at or above 
the knee and the method is solely based on the measurement of extensive air 
showers, i.e. the secondary particles created in the atmosphere and detected 
by large arrays on the ground. 

The origin of the knee has been one of the fundamental problems of cosmic-ray 
physics, and it has been discussed for decades. Several models have been pre¬ 
sented predicting different composition at the knee energies, and could only be 
identified by the experimental evidence on the composition. Some new exper¬ 
imental efforts have been devoted to the study of cosmic rays in recent years. 
These experiments are based, for example, on multi-parameter measurement 
of extensive air-showers, on shower maximum measurement by Cerenkov or 
fluorescence detectors and on underground multimuon measurements (see, for 
example, Ref. [1] and references therein). Their conclusions, however, have so 
far been diverse, implying the need for further studies, especially using differ¬ 
ent approaches. The results are also known to be strongly model dependent. 

EMMA (Experiment with MultiMuon Array) uses a different approach. It 
is not the first underground cosmic-ray experiment (see, for example, Refs. 
[2,3,4,5,6]), but it differs significantly from previous underground experiments 
with its ability to measure the lateral distribution function of underground 
muons. In EMMA the composition analysis is based on the lateral distribu¬ 
tion of high-energy muons and on their multiplicity. The muons detected by 
EMMA are generated in the upper part of the air shower close to the primary 
interaction. 


2 Experimental details 


The present experiment will be carried out at the depth of 85 metres in the 
Pyhasalmi mine (owned by the Inmet Mining Corporation Ltd., Canada), 
which is situated in the middle of Finland. 
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Fig. 1. Schematic layout of the EMMA array in the underground site at 85 metres. 
The array consists of nine units (red squares) each having an area of about 15 m 2 . 
The horizontal and vertical dimensions are approximately 50 and 20 metres, respec¬ 
tively. The three crossed units occupy detectors in two layers for shower direction 
measurement. The gaps between the units can be filled with detectors and the two 
open squares indicate the possible additional units if the array is expanded. 

The detector array is placed at the depth of 85 metres (corresponding 240 
m.w.e) which gives a threshold energy of muons of approximately 50 GeV. 
The rock overburden filters out all other particles of the air shower except the 
high-energy muons (and neutrinos). A schematic layout of the array is shown 
in Fig. 1, where each of the nine units have an area of about 15 m 2 . The area 
between the units is close to the unit area. 

The EMMA detector array consists of drift chambers (muon barrel (MUB) 
detectors) previously used in the LEP-DELPH1 experiments at CERN [7]. 
The setup is able to measure the arrival direction of the air shower, the muon 
multiplicity and their lateral distribution. This is carried out by an array of 
separate detectors with the total detector area of about 150 m 2 . Part of the 
array is in two layers with a vertical distance of 2.5 metres between the layers 
in order to obtain the direction information (see Fig. 1 for details). 

Most of the drift chambers have an active volume of 365x20x1.6 cm 3 . A 
plank consists of seven chambers (partly overlapped) having an area of ap¬ 
proximately 3 m 2 . The layout of Fig. 1 requires 60 planks. The drift chambers 
operate in the proportional mode, with Ar:CC >2 (92:8) nonflammable gas mix¬ 
ture. Due to the safety issues, the gas mixture does not contain CH 4 . Each 
drift chamber can provide up to three signals, one anode signal and two delay 
line signals (near and far), which can be used to localise the positions of par¬ 
ticles passing through the chambers. The measured position resolution of the 
chambers in our setup is about 1 cm and 3 cm along the drift direction and 
along the drift line, respectively. 

Recycling the old muon barrel detectors provides a the possibility to built the 
array at low costs. Also the use of existing caverns of the mine reduces the 
costs. The data-acquisition electronics is build by ourselves or bought as new. 
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Fig. 2. Eight muons hitting two planes of detectors with randomly chosen shower 
direction and muon hit positions. The lines connecting the hits in the upper and 
lower planes result from the tracking code. 

The underground array is under construction and the data recording can be 
started around middle of 2006 with a partial-size array. The full-size array is 
expected to be in operation by the end of 2007. 


3 Data analysis 

3.1 Muon tracking 


An important part of the data analysis software is the muon tracking routine. 
As it is obvious that an observed shape of the shower partially depends on the 
angles it hits the detectors, the shower can only be reconstructed by determin¬ 
ing the angles of individual muons. In EMMA the shower reconstruction will 
be carried out by tracking muons as they hit simultaneously to two detector 
planes (the three double-layer units shown in Fig. 1). 

The simulated performance of the muon tracking routine is illustrated in Fig. 2 
using an example where eight simultaneous muons penetrate two detector 
planes with a randomly chosen shower angle. The tracking code employs an 
automatic routine which connects the pairs in two planes and rejects pairs too 
far from the average angle. 

The tracking routine will be tested with real cosmic-ray muon data at the 
surface in the end of 2005 and in the beginning of 2006 using system similar 
to Fig. 2, and it will be further developed. 
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Fig. 3. Simulated (CORSIKA & QGSJET [8,9]) lateral distributions of muons for 
1, 3 and 10 PeV proton- and iron-initiated showers. The cutoff in muon energy is 
50 GeV. 

3.2 Locating core position 


In order to reconstruct the lateral distribution of underground muons and 
their multiplicity, the core of the air shower has to be located accurately (i.e. 
within metres). A two-dimensional fit routine was developed for the extrac¬ 
tion of the core location. The routine was tested with simulated air-shower 
samples. Depending on the hit position of the core, it can be extracted with 
an accuracy better than three metres in most cases if the core hits the detec¬ 
tor unit or between them. If the core is elsewere, the accuracy is weaker and 
the distribution slightly biased, but this can be corrected for. In the current 
version of the core-location software the routine uses model-dependent lateral 
distributions. However, we are developing a routine that works in the model- 
independent way, which would allow a direct comparison between measured 
lateral distribution and different model predictions. 


3.3 Extracting composition 


The multiplicity measurements at distances between 20 and 30 metres from 
the core position are expected to be the optimal cases for separating different 
cosmic-ray primary particles. As the core is located and the lateral distribution 
is fitted for each shower, the measured lateral distributions are used to extract 
the composition information as a function of the primary energy. The number 
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Muon density at 1 m 

Fig. 4. The separation of primary cosmic-ray composition and energy, in 
two-dimensional plot of p{ 1 nr) versus p (20 m). See text for details. The plot is 
based on simulations of proton and iron primaries with energies of 1, 3 and 10 PeV. 

of muons at the core can be used as an indicator for the primary energy. The 
energy resolution of EMMA according to this preliminary analysis is somewhat 
moderate. 


4 Expected results 


As no experimental results have been obtained yet, expected results based on 
simulations are introduced in this section to present some of the performance 
of the EMMA array. The simulations serve, at the same time, also a way to 
develop and further improve the methods to analyse and to interprete the 
data. The air-shower simulations and particularly the data analysis methods 
are still preliminary but they already indicate the good capabilities of the 
EMMA array. 

According to cosmic-ray air-shower models (e.g. CORSIKA & QGSJET [8,9]) 
the lateral distribution, or its gradient, of high-energy muons is sensitive to 
primary cosmic-ray particles and their energies, and to hadronic interactions. 
Also the muon multiplicity is sensitive to the energy. The simulated (using 
CORSIKA & QGSJET) lateral distribution functions of muons are shown in 
Fig. 3 for three discrete energies of 1, 3 and 10 PeV for proton- and iron- 
initiated showers, being the most relevant for the present experiment. Only 
muons with energies above 50 GeV are included in Fig. 3. 
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Primary energy [PeV] 


Fig. 5. Integrated shower rates above given primary energies for two effective de¬ 
tection areas. 


The method of the cosmic-ray composition study in the present work is illus¬ 
trated in Fig. 4 where the separation of two primary cosmic-ray particles with 
three different energies is shown. The figure is generated by fitting the lateral 
distribution functions to muons in an air shower and extracting the muon den¬ 
sities at distances of 1 and 20 metres (p( 1 m) and p(20 m), respectively) away 
from the located core position. The composition and energy of the primary 
cosmic ray can then be extracted, for example, using variables p(l m), which 
is related to the energy, and p(20 m), which is related to the composition. 


In reality, different groups (or groups of elements) in the plot would be less 
pronounced than those simulated and shown in Fig. 4. This is due to poorer 
statistics and an uncertainty related to the core position determination, and 
also continuous cosmic-ray primary energy spectrum (Fig. 4 shows only three 
discrete energies for only two primary particles). 


The effective area of the array determines recorded statistics. The area in the 
EMMA experiment is expected to be between 100 m 2 and 1000 m 2 . Fig. 5 
shows the expected numbers of air showers collected as a function of cosmic- 
ray primary energy. 
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5 Summary 


A new underground cosmic ray experiment EMMA is under construction and 
it is expected to start recording data in the full scale by the end of 2007. 
With a partial-size array the data recording can be started by the middle of 
2006. The analysis of simulated air showers shows that the primary cosmic-ray 
composition could be resolved (with a two-component model) at and above 
the knee energies. A possibility for the model-independent way of the deter¬ 
mination of the muon lateral distribution would allow to improve high-energy 
interaction models. Due to new method used, the EMMA experiment could 
provide comprehensive (and perhaps new) information on the composition of 
the cosmic rays at the knee region within the next few years. 
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